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ABSTRACT 

This paper proposes a class of estimators for population correlation coefficient 
when information about the population mean and population variance of one of the 
variables is not available but information about these parameters of another variable 
(auxiliary) is available, in two phase sampling and analyzes its properties. Optimum 
estimator in the class is identified with its variance formula. The estimators of the class 
involve unknown constants whose optimum values depend on unknown population 
parameters.Following (Singh, 1982) and (Srivastava and Jhajj, 1983), it has been shown 
that when these population parameters are replaced by their consistent estimates the 
resulting class of estimators has the same asymptotic variance as that of optimum 
estimator. An empirical study is carried out to demonstrate the performance of the 
constructed estimators. 
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1. Introduction 

Consider a finite population U= {1,2,..,1,..N}. Let y and x be the study and auxiliary 
variables taking values y; and x; respectively for the ith unit. The correlation coefficient 
between y and x is defined by 


Py = Syx I(SySx) (1.1) 
where 
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Based on a simple random sample of size n drawn without replacement, 
(xi, yi), 1 = 1,2,...,n; the usual estimator of p,, is the corresponding sample correlation 


Cea 


Ms: 


I 
un 


i 








coefficient : 
T= Syx/(SxSy) (1.2) 
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where s =(n-1) ¥ (y; - Xx; -¥), s2 =(2-1) (x, - 2) 


i=l i=l 
s? =(n-1)" Èo- y) y= mY k= n'a. 


The problem of estimating p,, hes been EA taken up by various authors including 


(Koop, 1970), (Gupta et. al., 1978, 1979), (Wakimoto, 1971), (Gupta and Singh, 1989), 
(Rana, 1989) and (Singh et. al., 1996) in different situations. (Srivastava and Jhajj, 1986) 
have further considered the problem of estimating ,, in the situations where the 


information on auxiliary variable x for all units in the population is available. In such 
situations, they have suggested a class of estimators for p „ which utilizes the known 
values of the population mean X and the population variance S? of the auxiliary variable 


x. 
In this paper, using two — phase sampling mechanism, a class of estimators for 


Px in the presence of the available knowledge (Z and S2) on second auxiliary variable z 


is considered, when the population mean X and population variance S? of the main 
auxiliary variable x are not known. 


2. The Suggested Class of Estimators 

In many situations of practical importance, it may happen that no information is 
available on the population mean X and population variance S? , we seek to estimate the 
population correlation coefficient »,, from a sample ‘s’ obtained through a two-phase 
selection. Allowing simple random sampling without replacement scheme in each phase, 
the two- phase sampling scheme will be as follows: 
(i) The first phase sample s* (s° cU ) of fixed size n,, is drawn to observe only x in 


order to furnish a good estimates of X and S2. 
i) Given s*ř, the second- phase sample s (scs*) of fixed size n is drawn to 


observe y only. 
Let 


= (Yn) dix, = (Yn) Z y: x =m) Ym s =a- Z-a), 


ies ies ies* ies 


s? =(n, -1'5 (x; a ; 


We write u=x/x* ,v =s? / s? . Whatever be the sample chosen let (u,v) assume values in 
a bounded closed convex subset, R, of the two-dimensional real space containing the 
point (1,1). Let h (u, v) be a function of u and v such that 
h(1,1)=1 (2.1) 

and such that it satisfies the following conditions: 
1. The function h (u,v) is continuous and bounded in R. 
2. The first and second partial derivatives of h(u,v) exist and are continuous and 

bounded in R. 

Now one may consider the class of estimators of p,, defined by 

Pra =r hlu, v) (2.2) 

which is double sampling version of the class of estimators 


7, =r f (u*,v*) 


Suggested by (Srivastava and Jhajj, 1986), where u* =x/X, v* =s? /S? and (x, 5?) are 
known. 

Sometimes even if the population mean X and population variance S? of x are 
not known, information on a cheaply ascertainable variable z, closely related to x but 
compared to x remotely related to y, is available on all units of the population. This type 
of situation has been briefly discussed by, among others, (Chand, 1975), (Kiregyera, 
1980, 1984). 

Following (Chand, 1975) one may define a chain ratio- type estimator for p,, as 
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where the population mean Z and population variance $? of second auxiliary variable z 
are known, and 


z = (I) Xi > s? =(n, Ek az 
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are the sample mean and sample variance of z based on preliminary large sample s* of 
size n;(>n). 


The estimator J; in (2.3) may be generalized as 
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where a@,'s (i=1,2,3,4) are suitably chosen constants. 





Many other generalization of „is possible. We have, therefore, considered a 
more general class of p, from which a number of estimators can be generated. 

The proposed generalized estimators for population correlation coefficient p,, is 
defined by 


Pig =Vt(U,V, W, a) (2.5) 
where w=7"/Z ,a=s*/S? and t(u,v,w,a/ is a function of (u, v, w,a) such that 

t (1,1,1,1)=1 (2.6) 
Satisfying the following conditions: 
(i) Whatever be the samples (s and s) chosen, let (u,v,w,a) assume values in a closed 
convex subset S, of the four dimensional real space containing the point P=(1,1,1,1). 
(11) In S, the function ¢(u,v,w,a) is continuous and bounded. 
(iii) The first and second order partial derivatives of t(u,v,w, a) exist and are 

continuous and bounded in S 

To find the bias and variance of 6,, we write 
sy =S)(+e,),%¥=X(+e,),x%° =X +e) si =S] (+e) 
sy =Si(l+e3),2 =Z +e), s? =S7(1+e4),8,, =Sy(+e,) 
such that E(eo) =E (e))=E(e2)=E(es)=0 and E(eř)=0 Y i= 1,2,3,4, 
and ignoring the finite population correction terms, we write to the first degree of 
approximation 


E e? )= (S -1)/n, Ele? )= C? /n, Ele}? )= C2 /n, Ele? )=(Sou -1)/n, 
rle? )=(8 640 —1)/n, E(e3?)= C? fn .E(e??)= (Sos —1)/n, > 

E(e 2)= fe wn /n, E( E(e, €,)=5yC, /n, Eley et )=5y9C,/m,, 
E( 


eo €>) =(Sy9 —L)/n, Ele, e} )= (Sz —1)/n,, lezez) Dai Csi Nis 
Ele, ei)=(8 (S39) =1)/n, -E Eley es)= f8s10/2x)- 14/7, 
Ele, e*)= C? Jn, Ele e,)=8 S0C, /n, Ele, ez)= 5C, /n,, 
Ele, e} )= PEC Cnet Ele, e3)=dynC,/n, Ele, €5)=(By9C./ Py )/ns 
Ele; e,)= Ome im Ele e; )= 30C /n Ele e3)= p,.C,C. Jn, 
Ele; ej)= SoC, /n, Ele} e Je (5,20C./Pyx) iis 
Ele, e3)= (Sou -1)/n, Ele, eż)= SonC.[n, Ele, e}) (Sox. -1)/n,, 
E(e, e; )= (6 a] i\/n, Ele; e} eż )= 5m1, /n,,; 
Ele; e7)= (ĉn -1 1)/n, Ele} e s= (60/2) -1/7 
Ele; e*)= wC; /n E Ele; €s)=(Cz/Py)/m, 
Ele; e;)= Q 12/2 )- I}/n,. 


where 
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non-negative integers. 

To find the expectation and variance of /,,, we expand t(u,v,w,a) about the point 
P= (1,1,1,1) in a second- order Taylor’s series, express this value and the value of r in 
terms of e’s . Expanding in powers of e’s and retaining terms up to second power, we 
have 


E( pu =P +o(n") (2.7) 


which shows that the bias of 4,, is of the order nand so up to order n’ , mean square 
error and the variance of ,, are same. 


Expanding (Da =p l; retaining terms up to second power in e’s, taking 
expectation and using the above expected values, we obtain the variance of 6,, to the 
first degree of approximation, as 
Var( Ôa) = Var(r) + (py, (MLC ity (P) + (Soa — Dt (P) — At, (P) — Bt, (P) + 2803C ,t (P)ty (P)] 

=(P INDICE (P) + (Soa — Dt (P) - CZt3 (P) — (o4 — Vtg (P) - At, (P) - 
Bt, (P) + Dt, (P) + F t, (P) + 25 q39C ,t (P)ty (P) — 2803C t; (P)t,(P)] 
(2.8) 
where 1¢7(P), t(P), t3(P)and t4(P) respectively denote the first partial derivatives of 
t(u,v,w,a) white respect to u,v,w and a respectively at the point P= (1,1,1,1), 


Var(r)= AA IME | p?,)+ (1/ 4)(Sq49 + F400 + 26.299) — {(6139 + 9310) / Pyxt] 9) 
A= {0419 + 6039 — 200120 / P yx) IC x» B = {6209 + Bou — 2(6130 / Py) I> 
D = {őz + 6021 — 206111 / Pyx IC, F = {8x02 + O22 — 206112 / Py) } 


Any parametric function t(u,v,w,a) satisfying (2.6) and the conditions (1) and (2) can 
generate an estimator of the class(2.5). 


The variance of 6,, at (2.6) is minimized for 
[A(S = 1) zz BO ox C, ] 
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Thus the resulting (minimum) variance of 6,, is given by 
2 A -B 2 
min.Var(p,,)=Var(r)— (2-4 ohl 2 -+ a [Ca Poo } 
t- Mi 4C; 4(S 40 — Soo 7 D (2.11) 
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It is observed from (2.11) that if optimum values of the parameters given by 
(2.10) are used, the variance of the estimator ,, is always less than that of r as the last 
two terms on the right hand sides of (2.11) are non-negative. 

Two simple functions t(u,v,w,a) satisfying the required conditions are 


t(u,v,w,a)= 1+a,(u-1)+a,(v-l)+a,(w-1)+a,(a-l 
t(u,v,w,a) =u" vwa" 
and for both these functions t,(P) =a@,, t2 (P) =a@,, t3 (P) =a@, and ty (P) =a@,. Thus one 
should use optimum values of @,,@,, @,and a,in p,, to get the minimum variance. It is 
to be noted that the estimated 6,, attained the minimum variance only when the optimum 
values of the constants œ, (i=1,2,3,4), which are functions of unknown population 
parameters, are known. To use such estimators in practice, one has to use some guessed 
values of population parameters obtained either through past experience or through a 
pilot sample survey. It may be further noted that even if the values of the constants used 
in the estimator are not exactly equal to their optimum values as given by (2.8) but are 
close enough, the resulting estimator will be better than the conventional estimator, as has 
been illustrated by (Das and Tripathi, 1978, Sec.3). 

If no information on second auxiliary variable z is used, then the estimator ô, 
reduces to 6,, defined in (2.2). Taking z = 1 in (2.8), we get the variance of 6,, to the 


first degree of approximation, as 
Var D,a) =Vanr) + P lc2n? (1,1) + (S49 — DAZ 1) — Ah (1,1) — Bh, (1,1) + 2603C h Dh, an] 


(2.12) 
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. bl) = ABC, AO eal) (2.13) 
2C, (Soso — 6430 -1) 


hı(1,1) = 
Thus the minimum variance of 6,,, is given by 
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It follows from (2.11) and (2.14) that 
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which is always positive. Thus the proposed estimator ô, is always better than ,,- 


min. Var( Ôu )-min.Var( yy )= (2 /m) l ] (2.15) 


3. A Wider Class of Estimators 
In this section we consider a class of estimators of »,, wider than ( 2.5) given by 


Pea =8(r, u, v,w,a) (3.1) 


where g(r,u,v,w,a) is a function of 7,u,v, w,a and such that 


g(p,1,1,1,1)= >, and [22] =1 
(p,1,1,1) 


or 





Proceeding as in section 2, it can easily be shown, to the first order of approximation, that 
the minimum variance of Pea is same as that of 9,, given in (2.11). 

It is to be noted that the difference-type estimator 
ra= r + a, (u-l) +a, (v-1) + a, (w-l) + æ, (a-1), is a particular case of Pea , but it is 


not the member of /,, in (2.5). 


4. Optimum Values and Their Estimates 


The optimum values t)(P) = a, t(P) = $ ,t{P)=y and t{P)=8 given at 
(2.10) involves unknown population parameters. When these optimum values are 
substituted in (2.5) , it no longer remains an estimator since it involves unknown 
(a,f8,yv,6), which are functions of unknown population parameters, say,, ô pym (P, qM= 


0,1,2,3,4), C» C, and p,, itself. Hence it is advisable to replace them by their consistent 


estimates from sample values. Let (a, B, 7,6 ) be consistent estimators of t)(P),t.(P), 
t;(P) and t4(P) respectively, where 
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(4.1) 
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= Sfi + Soro — 26199 (TIC, » B =[ő + Oo — 2(O139/ r], 
5 Soy + dpi —2(6,,/ r)IC F = [S22 + o2 — UÊ /r)], 
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We then replace (a, 8,7 ,8 ) by (â, 8,?,ô ) in the optimum ô, resulting in the estimator 
Ê say, which is defined by 


A A 


Py =rt (u,v,w,a,Ĝ „B®, (4.2) 


where the function t*(U), U= (u,v,w,a,@, Ê, 7,6) is derived from the the function 


t(u,v,w,a) given at (2.5) by replacing the unknown constants involved in it by the 
consistent estimates of optimum values. The condition (2.6) will then imply that 























r(P*)= 1 (4.3) 
where P*= (1,1,1,1, @,8,7,6) 
We further assume that 
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Expanding t*(U) about P*= (1,1,1,1, a, 2 ,y ,ô ), in Taylor’s series, we have 


Py =r[t (P*)+ (u -Dt (P) +0 -DE (P") +(w- D (P*)+ (a- ty (P*)+ (Â -ati (P’) 
+ (B — Pti (P*) + (? — yt (P*)+ (6 — Ôt; (P*) + second order terms] 
(4.5) 
Using (4.4) in (4.5) we have 


Py =rlt+(u-Dat+(v-DB+(w-Dy + (a —1)8 + second order terms] (4.6) 


Expressing (4.6) in term of e’s squaring and retaining terms of e’s up to second degree, 
we have 


A 1 * * * * 
(Pa -p = ppl es -ey -e,)+aæle —e,) + Ble, -e3) tye; +e] (4.7) 


Taking expectation of both sides in (4.7), we get the variance of p} to the first degree of 
approximation, as 





A? {CAIC So ~ BY’ 
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which is same as (2.11), we thus have established the following result. 


Var( 4) =Var(r) c De 
1 (4.8) 
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Result 4.1: If optimum values of constants in (2.10) are replaced by their consistent 
estimators and conditions (4.3) and (4.4) hold good, the resulting estimator ô% has the 
same variance to the first degree of approximation, as that of optimum ô, . 


Remark 4.1: It may be easily examined that some special cases: 


a oat O â Â PÂ ee o {1+ &@u-1)+7(w-D} 
(i) Pidi ru v w f (ii) Praz “A Bw 1) bla 1} 





(iii) pis =r [1+ @u-l) + Bu-1) + P(w-1) + d(a-D)] 





(iv) Puy =r- âu -1)- Buu-1) -P(w-1) -d(a-D]" 
of ô satisfy the conditions (4.3) and (4.4) and attain the variance (4.8). 


Remark 4.2: The efficiencies of the estimators discussed in this paper can be compared 
for fixed cost, following the procedure given in (Sukhatme et. al., 1984). 
5. Empirical Study 

To illustrate the performance of various estimators of population 
correlation coefficient, we consider the data given in (Murthy, 1967, p. 226]. The 
variates are: 
y=output, x=Number of Workers, z =Fixed Capital 
N=80, n=10, n;-25, 
X = 283.875, Y =5182.638, Z =1126, C, =0.9430, C, = 0.3520, C, = 0.7460, 
O93 = 1.030, ów, = 2.8664, pı =1.1859, p =3.1522, o =1.295, mo = 3.65, 
Oo) = 0.7491, Ôi = 0.9145, ô = 0.8234, O49 = 2.8525, 
Ôn = 2.5454, z =0.5475, 549 = 2.3377, 645, = 0.4546, p = 2.2208, z% = 0.1301, 
O49) = 2.2667, p,, = 0.9136, Px, = 0.9859, p = 0.9413. 
The percent relative efficiencies (PREs) of ĝia» Ênras Êa With respect to conventional 
estimator r have been computed and compiled in Table 5.1. 


Table 5.1: The PRE’s of different estimators of p 





Estimator r Pha Pag (OL Puy) 




















PRE(..r) 100 129.147 305.441 





Table 5.1 clearly shows that the proposed estimator „(or ô% )is more efficient 
than rand /,,. 
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